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We study the scaling property of the ratio of the critical temperature T c to the square root of the string tension 
a in the SU(3) pure gauge theory using a renormalization group (RG) improved action. We first determine the 
critical coupling j3 c on lattices with temporal extension Nt = 4 and 6, and then calculate the static quark potential 
at the critical couplings on lattices at zero temperature. The values of T c /^/a in the infinite volume limit are 
identical within errors, while they are slightly larger than the value extrapolated to the continuum limit with the 
standard action. We also note that the rotational invariance of the static quark potential is remarkably restored 
in the both cases, and that the potential V(R) in physical units scales in the whole region of R investigated. 



1. Introduction 

In studies of lattice QCD, it is important to 
control and reduce finite lattice spacing effects. 
Several improved actions have been proposed for 
this purpose and some of them have been tested 
for the scaling behavior of the critical temper- 
ature T C [|IJ. In this work we study the scaling 
property of T c /^/a using a RG improved action 
proposed by one of the authors 0|: 

S g = l/g 2 {c y^(plaquette) + ci ^(1 x 2 loop)} 

with c\ = —0.331 and cq = 1 — 8ci. 

The outline of the calculation of T c j ^fu is as 
follows. First we determine the critical coupling 
ft's on 12 3 x 4 and 18 3 x 6 lattices. Then the 
string tensions at the two /3 c 's are evaluated from 
smeared Wilson loops on 12 3 x 24 and 18 3 x 36 
lattices, respectively. The string tensions in the 
infinite volume limit are estimated for the both 
cases, with a finite size scaling analysis applied 
using the result of f3 c on a 16 3 x 4 lattice. 

2. Determination of f3 c 

The critical couplings /3 c 's for the finite temper- 
ature phase transition in the SU(3) pure gauge 
theory are determined as the location of the peak 

*Talk presented by T. Kaneko at Lattice 96. 



Table 1 

Parameters of the finite temperature simulations 



a 


sweep 


thermalization 


12 3 x 4 






2.250 


12000 


2000 


2.275 


125000 


40000 


2.300 


10000 


1500 


16 3 x 4 






2.283 


220000 


40000 


2.290 


240000 


40000 


18 3 x 6 






2.5000 


120000 


15000 


2.5125 


256000 


50000 


2.5250 


210000 


60000 


2.5375 


135000 


5000 



of the susceptibility of the Z(3) rotated Polyakov 
line. Parameters of the simulations are summa- 
rized in Table [|. Gauge fields are updated by 
Cabibbo-Marinari-Okawa pseudo-heat bath algo- 
rithm here and also for the calculation of the 
string tension. The results of the susceptibility 
with the histogram method applied on the 12 3 x 4 
lattice are shown in Fig. [|. The error of (3 C is de- 
termined by the jack knife method with bin size 
of 3000 sweeps for the 12 3 x 4 and 18 3 x 6 lattices 
and 6000 sweeps for the 16 3 x 4 lattice. We have 
checked that the errors are stable for the bin size 
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Figure 1. Susceptibility on the 12 3 x4 lattice. 
Filled symbols represent the location of simula- 
tion points. 

larger than these values. We obtain 



2.2827(16) 
f3 c = { 2.2863(10) 
2.5157(7) 



(12 3 x 4) 
(16 3 x 4) 
(18 3 x 6) 



(1) 



3. Evaluation of the string tension 

We evaluate the string tension on a 12 3 x 24 
(18 3 x 36) lattice at j3 = C for the 12 3 x 4 
(18 3 x 6) lattice. After at least 5000 ther- 
malization sweeps, we measure Wilson loops 
at every 200 sweeps. The spatial paths of the 
loops are taken to be multiples of spatial vectors 
(1,0,0), (1,1,0), (2,1,0), (1,1,1), (2,1,1), (2,2,1). 
The number of configurations is 200 (100) at (3 C 
for the 12 3 x 4 (18 3 x 6) lattice. We carry out the 
smearing technique proposed in Ref.Q up to 30 
(40) times on the 12 3 x 24 (18 3 x 36) lattice. The 
potential V(R) and the overlap function C(R) 
are extracted by fitting the Wilson loops to 



W(R, T) = C{R) exp [-V{R)-T] . 



(2) 



The optimized smearing time is determined from 
the condition C(R) ~ 1.0. The string tension is 
determined by fitting V(R) to 



V(R) = V -a/R + aR. 



(3) 
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Figure 2. Data for the potential and its fitting 
curve on the 12 3 x24 (above) and 18 3 x36 (below) 
lattices. 

We make correlated fits for the both fits. Fitting 
ranges are 



T = 3-5 i? = V6-4V2 (12 3 x24) 
T = 4 - 7 i? = 2V3 - 4\/5 (18 3 x36) 



(4) 



These fitting ranges are determined by \ 2 /df °f 
the fits and the stability of the fitted values of 
a and a. Note that the changes of the fitting 
ranges are consistent with the change of the scale 
between j3 = j3 c {N t = 4) and f3 c {N t = 6). 

The values for V(R) are plotted in Fig. |^, which 
are excellently fitted by the rotational invariant 
fitting curve, EqJ^. The averaged deviation is less 
than 0.5% on the both lattices. The results of Vq, 
a, a are summarized in Table 0. 



3 



Table 2 
Results of Vb , 



M 
4 
6 



V 

0.630(20) 
0.627(18) 



0.295(14) 
0.297(19) 



a 

0.1493(25) 
0.0655(12) 
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Figure 3. V/T c vs. i?-T c . The constant term in 
the potential is fixed so that both curves have the 
same value at R-T c = 0.866. 



4. Scaling properties 

In Fig. ||, the values of V(R)/T C are shown as a 
function of R-T c . We note that the data obtained 
on the both lattices are in excellent agreement in 
the whole R-T c region. 

We obtain the values of T c /y/a in the infinite 
volume limit: 



T, 



0.647(- 
0.651(- 



-7)(5)(3) 
-8)(6)(6) 



(N t = 4) 
(N t = 6) 



(5) 



The error in the first bracket is due to the shift of 
the critical coupling from (3 C at N s /N t = 3 to the 
one in the infinite volume limit the second 
one is the statistical error and the third one is 
due to the error in the values of 

The error in the first bracket above has been es- 
timated as follows: We have /3 C °° = 2.2889(29) and 
2.5219(52) for Nt = 4 and 6 lattices, respectively, 
assuming both for Nt = 4 and 6 lattices a finite 
size scaling relation [Q] obtained from the data of 
/3 C on the 12 3 x 4 and 16 3 x 4 lattices: (3 c (N t , oo) = 
Pc{N t ,V) + 0.168(121) • N t 3 /V. Then the values 
of the string tension at /3^° are estimated assum- 
ing an exponential scaling *J~o — Aexp(— B0). 
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Figure 4. T c /a x / 2 vs. 1/N t 2 . All symbols repre- 
sent the values in the infinite volume limit. 



Our results (||) are shown in Fig. |i| together 
with the results in Ref.|5) using the standard 
action. Our values are slightly larger than the 
value extrapolated to the continuum limit with 
the standard action 0.629(3). Using the aver- 
age value of (||) and the experimental value a = 
(420MeV) 2 , we estimate T c - 276(3) (2)MeV. We 
note that a -1 w 1.10 GeV and 1.66 GeV, respec- 
tively, in the Nt = 4 and 6 cases, when we use 
the value of T c in a" 1 = N t T c . 

Numerical simulations are performed with Fu- 
jitsu VPP500/30 and HITAC H6080-FP12 at the 
University of Tsukuba. This work is in part 
supported by the Grants-in-Aid of Ministry of 
Education, Science and Culture (Nos.07NP0401, 
07640375 and 07640376). 
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